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ABSTRACT 


The primary objective of this study was to provide 
a contribution to the phenomeron "Suppression" as an 


aspect within the military environaent. 


Analytical models explaining these aspects were 
developed in order to identify the influences to 
Suppressicn. Techniques are examined for including tne 
Suppressive effects of weapon systems in Lanchester 
type combat models, which may be useful in wargame 
evaluations of military judgements, and in force level 
planning. The study also provides techniques to 
analyze and fit experimental data to the analytical 


models. 


The data to verify the models were cbtained fron 
related experiments performed by Combat Development 
Experimentation Command (CDEC), Fort Ord, California. 


The result for the modelling approach to 
Suppressicn indicates source dependences on 


guantitative as well as on qualitative features. 
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A. INTRODUCTION 


In recent years considerations of Human Factors asrects 
in the military has gained more and more importance. The 
study of psychological, physiological, and environmental _ 
conditions and their influences on the performance of men in 
IMan-Machine-Systens" supports the develcpment of new 
doctrines, design cf weapon systems as well as training 


programs for troops. 


One phenomenon in military man-machine-systens is 
suppressicn, Modelling this phenomenon has only recently 


been given much attention. 


This is partly kecause modellers did not understand 
exactly neither the causes of suppression ner how it affects 


the course of combats. 


There is an intuitive feeling that when a soldier or a 
combat unit is being fired upon, it will be less effective 
than when it is not receiving tire. This is generally 
referred to as suppression, but it can include much more. In 
the broadest sense, suppression can effect individuals, 


units, cr weapon systems in different types cf combat. 


This paper will limit itself to individuals or small 
infantry units. 





——— — 


In order to be able to build up functional relationships 
based on the suppression idea, between different categories 
of people and time, this paper will first exanine different 
rates as input to analytical Lanchester-type-models. The 
rates will be based on stochastic aspects. The modelling is 
done under different viewpoints; their results will be 


ccnpared and discussed. 


In the second part of the thesis, the parameters cf the 
models will te tested in a regression analysis against real 
world data which were placed at the writers' disposal py 
C DEC, 


The goal of this paper is not tc present a final 
framework cf suppression; however it may contribute to 
Clarification of aspects of suppression, and help to ented 


it as a component in futurelarge size models. 
B. BACKGROUND 


The strains imposed on individuals in our society are 
constantly increasing. Modern technologies and constant 
efforts for improvements of standards of living lead to 
growing difficulties in adjustments or sometimes to complete 
failures to adjust. These facts create stress and we may 
cbserve that the degree of stress increases with the 


difficulty of the adjustment-problem. 


The term stress will be used as a substitute for what 
might be called otherwise as anxiety, conflict, emotional 
distress, extreme environmental conditions, ego-threat, 
frustration, threat to security, tensicn or probably 
arousal. Stress can be thought of as the result of almost 


any envircnmental interference [Appley,p.2). 








The stress-generating features of the civilian 
environment are great, but the environment created by modern 
warfare possesses additional features which result in an 
increase of stress. The combat environment created ty the 
weapcn power of tne enemy causes a constant threat to life. 
The soldier has to operate under this threat and naturally 
he will respond with ccnstantly recurring fear. This may 
break down the soldier's psychological and physiolcgical 


resistance. 


Fear and anxiety in battle is coamon, being experienced 
cy between 80 and 90 % of combatants. Pains in the stonmacnh, 
fatique, dizziness, perspiration, and enhanced heart-teat 
are some vegetative correlates to fear and anxiety. Of 
course the moment when the individual soldier reaches his 
breaking point varies and depends en individual 
psychological and pnysiological resistance and the severity 
of the battle. 


One interesting observation from Agrell was that 
auditory sensations convey the stress of pattle most 
strongly and most directly. The psychcicgical effect of 
weapons is directly related tc their sound level and the 


frequency with which the sound occurs (Agrell,p.215]. 


Each enemy grenade causes tne soldier to react 
constantly with fear. Stress and fear can have a 
significant sensory operating characteristic, e. g. the 
detection threshold and/or sensitivity may decrease as a 
result of stress [Weltmann, G, Christianson, R.A. and 
Egstrom, G.H., pP. 423-430]. 


C. MILITARY SUPPRESSICN IN CCHBAT-ENVIRONNBENTI 








When weapons are used in combat, there are two types of 
effects that they have. The first type of effect is physical 
damage or injury to the target and the seccnd type of effect 


is psychological. 


This second effect of a weapon has led to the term or 
concept "Combat Suppression". There exist a general telief 
that fire suppression is important, but the importance of 


suppression effect on combat outcomes as compared to the 


effect cf other areas such as firepcwer, mobility, 
intelligence, command/control has not been quantified 
adequately. 


1. Definitions 


a-—D A-— A— EM — — ee ee ee — 


SUEEression can be generally defined as 
the temporary degradation ın the quality 
Of performance of an individual Zcläier 
cr unit by an internal cr external 
stimulus. 


CDEC; 1277 


A more useful operational definiticn in terms of 
performance capability changes is provided by the "Repcrt of 
the Army Scientific Adviscry panel Ad Hoc Group on Fire 


Suppression". 


Fire suppression is a process which 
causes temporary changes in  rperformagzance 
capabilities of the Suppressee 
(Suppressed system) from those expected 
when runctioning in an environment he 
KnowS to p@ fasSive. These changes are 
caused oy signals from delivered fire or 
the threat of delivered tire, and the 

result from pehaviors that are intende 

to lessen risk to the suppressee. 


Ad Hoc Group, 1976 


This definition emphasizes that suppression is not a 
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single effect which can be measured totally on a single 
quantitive scale. Suppression effects are multidimensional 


and the "amount of suppressicn" varies among these 
dimensions (e.g. fire impact pcints, soldier's 
characterists, and reaction to the rire, his combat 


experience etc.). 

Suppressive fire in a combat environment can Suppress a 
number of combat activities; for example: firing, search for 
and observation oz targets, movenents of units or command 


and Control, 


2. Structure Of Fire Suppression Procass 


Fire suppression is a complicated process invclving 
Bany physical, environmental, physiological, behavioral, and 
operational variables. The important point to emphasize is 
that the Eehavior involved is in response tc stimuli that 
criginate both externally (combat environment) and 
Internally (personal background, training and experience) to 
the soldier suppressee. These aspects, however, are not 
included in this paper. 

The intensity and duraticn of suppression can not be 
predicted from a knowledge of the combaz environment alcne. 
It requires an analysis of the underlying motivational and 
cultural factors and of the context of the combat 
environment. 

the Fig 1 [Ad Hoc Grcup, p.36] shows a schematic 
descripticn of a process when suppressive fire is delivered 


aud its affect on the combat result. 
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Ihe process is discribed by individual £uncticns. 


a. Signal Process 


The first process in fire suppression 1s the generaticn of 


Signals provided by suppressing weapons. 


Inputs to this process are the 


Characteristics of the weapon systems (caliker, 


amount of porpellant, warheadtypes etc) and 


the environmental characteristics (trajectory, 


platform, arrival points ecc). 


Characteristics of the weapon systems can vary in order to 


increase or decrease suppressive effects. 


Some 


parameters that are considered to be important to 


suppressicn signals are: 


Muzzle velocity (an increase in muzzle velocity is 


associated with an increase in signal variables); 


Caliber (as caliber increases, the firing signais 
and projectile Signals increase along with 
lethality); 


Erojectile weight ( penetration depends on weight 
and velocity at impact and increases shock cougling 


to ground) ; 


Warhead charge weight ( the explosive charge weight 


determines the energy in the pressure pulse); 


Additional parameters like fire frequency and 


proximity of shots could also be mentioned. 


Environmental characteristics are also variables but they 


can not re determined completely. 


Environment has an influence on the signal 
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generation -and transmission-process. For example, auditory 
Signals that result from the impact of prcjectiles depend 
heavily on the nature cf the object or material impacted. 

A soft yielding material such as dusty ground or 
sand receiving the impact of a projectile will produce a 
different pulse and sound than will hard unyielding ground 
under the sare impact. 

Sound signals can be attenuated Ey the shadowing 
effect or large obstacles or may be increased by echo or 
reverberation. Visual signals are strongly modified ty the 
condition of Lighting. Haze, fog, rain, and snow act 
Similarly to smoke and dust. Ihe visual field is also 
reduced and interrupted by terrain and other obstacles. 

The envircnment modifies the produced stimuli 
when they are transmitted to the location of the suppressee. 
Cutputs orf this process are the attenuated sensory signals 


that become imputs to tne hugan sensory receptors. 


Bs Human Process 


Many of tne determinants of the soldier's performance on 
the battle field are unkown or at least uncertain - thought 
of as influenced py chance factors. This emphasizes the 
difficulties of predicting human behavior in a combat 
environment. Tae human process (sensory and perception) 
converts the received signals into a perception of the risk. 

Eattle field stimuli effecting the individual are 
detected and converted into sensory data Ey a process such 
as vision and audition, so the sensory rrocess suggests that 
th2 weapon systems stimuli relevant to suppression are the 

~loudness, and 


-visual impact. 


There exist moderating factors that influence 
the operating characteristcs of the sensory process and that 


determines which stimuli are effective. 
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Sensor y modifiers (1:604 earplugs, 
night-vision-devices) serve to change users! sensitivity 
range. A major effect of these devices is to change the 
salience of stimuli. 

High concentration on an activity or a high 
level of effort on an activity (e.g. missile-gunner is 
tracking a target or reloading his system) may increase the 
absolute threshold. | 

The posture of a soldier (standing or sitting) 
and the sequence of posture (cbserving, ducking, observing) 
influences the sensory capabilities (e.g. observing for 10 
seconds continuously is not equivalent to observing 5 
seconds , ducking 10 seconds and observing 5 seconds). 

The perception process integrates sensory and 
other information into a perception or the risk. Risk refers 
tc the uncertainty of damage, injury, or loss. It 
characterizes decision situations in which the consequences 
of choosing an action are uncertain. 

If there is no uncertainty in the possible 
cutcome, there is no risk. 

Ferceived risk is a function of uncertainty and 
the subjective value the individual associates with each 
outcome. 

Ferceived risk represents the output of the 
combined ‘sensory and perception process. It depends cn tne 
individual's experience and training in assessing risk from 
sensory information. Also cover provided by the environment 


and the individual's pcsture may influence risk perception. 


C. Reaction Process 


Given the input perception of risk, this process causes 
physical and mental reactions, which depend on the 

“current mission 

-task 

-activity 
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-combat training doctrine and experience 
“group dynamics and 
-the quality of leadership. 


It is conjectured that two individuals who perceive the same 
degree of high risk , but who have different amounts of 
combat engagement experience, might be likely to react 
differently tc the risks. 

The soldier's reaction is also influenced ty his 
current state. A soldier who has recently ducked may be nore 
likely to duck than one who has not , given the same 
delivered fire. 

Prior reaction or seguence of reactions may te a 


good predictor of the coming reaction. 


d. Ferformance Effects Process 


Given the reactions of the human benavior process, it is 
conjectured that these directly afiect the ¡performance of 
certain activities of the surrressee in a calculable way. 

If for example the suppressee taxes cover, ne 
may fire less often and less accurately and also might be 
less vulneratle. The magnitude and duration of these changes 
in performance are dependent on the characteristics of the 
System employed by the suppresseé and the target of his 
detivity. 

So the nature and duration c£ change En 
performance capabilities is determined by the performance 


eftects process. 


3. Suppression In Field Experimentations 


— — — — — — — — — 


A fire suppression research progran requires 
significant experimentation on behavioral attitudes and 
reactions tc risk. This necessity causes *remendous 


difficulties in trying to induce actual behavior in soldiers 
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in field experiments. Former studies shcwed, that the 
soldiers felt true psychological stress only in situations 
in which they believed that they were in real danger. Such 
situations are difficult to contrive and te control. Social 
and ethical limits and legal regulations preclude tne 
introduction of actual physical risk. Scldiers must be 
taught the "rules and risk" defined in that context. The 
success Df playing the role, being an individual 
rarticipating in a combat engagement, depends on the 
scldiers' motivation and willingness. 

Because of these reasons and the  zultidimensional 
shape of the fire supfressicn process as mentioned earlier, 
Suppression in field experimentations may te restricted only 


to some variables involved in this process. 


D. APPROPRIATE OBJECTIVES 


The overall objective of a fire suppression research 
should be to relate changes in performance capabilities 
caused ky fire suppression [Ad Hoc Grcur, Ds 110]. 
kesponding more directly the following cbjectives may be 


determined: 


Indicating the effects of suppression on combat 
results, 1. e. to develop rates cf suppression. 
These values may be compared to other effected areas 
and probably employed in computer sirulations. These 
numbers represent two Kinds of variables: Weapon 
system variables and human suppressicn performance, 


given operational and environmentai conditions. 


Determining characteristics c£ suppressive fire 
Systems, characteristics which should be assigned to 
such a weapon system. Results are developed 
experimentally . Chapter IILI B. will display some 


evaluated parameters and constants tcr the developed 
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model, which may also represent the suppressive 


Characteristics of the used weapon systems. 


Reducing suppressive effects. Ways to reduce the 
effect of suppressive fire may also be considered as 
an appropriate objective of suppression research. 
Special training or equipment can be assigned to tae 
soldiers or new tactics can be developed. This 


cbjective is beyond of the research of this paper. 


ES POSSIBLz ALTERNATIVES 


In order to get information about the fire suppression 
process, previous investigations were based on interviews 
and questionaires, because valuable information of the fire 
suppressicn process is stored in the minds of combat 
veterans. Studies on veterans of the Vietnam conflict and 
the wars cf the Near zast would be especially useful, since 
newer weapons were employed and the combats were shorter and 
vcre intersive. 

These studies may provide a good insight tc the 
suppressicn process and/or may also deliver valuable inputs 


for the mcdelling approach. 
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À. ASSUMPTIONS AND SIMPLIFICATIONS 


Based on the foregcing discussion, the fcllowing model 
is a detailed model [I1Iaylor,1978,p.12] which starts cut by 
considering the behavior aspect of a numan being under tae 
influence of artillery fire power. It is assumed that such 
Simple ucdels, which represent a small part of a total 
scenario, can be used profitably to investigate systen 
dynamics of tore complex models. The value of the mcdel 
derives from the fact that it forms intuitively plausible 
and transparent subsets in a large compcsition of other 
subsets which determine the basic structure of the conplex 
cperational model. In other words, tne wnole is described 
in terms of the sum of its parts. 

The basic concern of the analytical model developed here 
“ill be to model the Eehavior of an individual experiencing 
artillery fire, considered as a function of time, where that 
behavior depends upon ammunition types and location of 
detonating rounds. 

On the basis or particular assumptions and 
Simplificaticns it will be possible to apply the results 
cEtained to a group of people (a force) cn a battelfield. 
The result of these considerations “ill provide a 
relationship between time and tne actual number of people 
affected ty the fire power of the artillery . This last step 
of the model is carried out by using Lanchester type 
equations, so called after the pioneering werk of F. q. 
Lanchester. Finally, the models enable one to estimate the 


total firepower of the force at any point in time; depletion 
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GE total firepower is caused by attrition and by 


suppressicn. 


Lo BATES FOR THE MODEL 


1. Rates Of Suppression 


The tasic considerations in the preceeding paragraph 
support the assumpticns, that the behavior of a suppressee 
can be expressed by a conditional probability of suppression 
as a function of miss distance r and aspect angle 8. This 
function is represented by the family of suriaces of the 


form: 


Pope expl- x UNE E 1 (2.1) 


where P(S/8,r) is the conditional probability of suppression 
given that a particular round impacts under a certain aspect 
angle © and a certain miss distance r away ircm the foxhole. 

The line along the angle 9=0° is identical with the 
line of sight. It is the main direction of ckservation. 

The constant K and “S are parameters, which are 
determined by the experiment itself and by the environmental 
conditions. 

Ihey can be influenced by £actors as discussed in 


chapter I which may be recalled here briefly. 


-type of ammunition 
-frequency of arrival of rounds 
-ferceptual damage 


-total time spent in the foxnole (learning process) 
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-noise appearance of the rounds 

-flash light intensity of the rounds 

-rerfornance 

-personal factors like age, personal condition, etc. 
-degree of stress 


-motivation 


[The computational evaluation of the constant K and e will be 
performed in chapter III B. and C. In particular it will be 
important to determine K aS it varies with different types 


erzsaumuniticn, 


Ihe mathematical conditions for the two parameters K and e 
are: 

0 <£ < 1 

K > O 

If 8 is held fixed ,0<9<360% and P(S/0,r) varies 
between O<P(S/6,r)<1 we obtain a family of functions, which 
is two-dimensional and shows an exponential relationship 
Estween F(S/98,r) and r. 

This is illustrated in the following figure, where 
tae angle $ is held fixed at C’and 180% The parameters K and 


£ are assumed as being 1500 and 0.7 respectively. 
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r/m = distance/meter 


P(S/60, ,r) 


1.0 

0.5 
a 75 7$ 50 100 — A 
r/m 2 r/m 
for 0, {80° foxhole for 6, 0° 


1 2 
Figure 2 - FUNCTION P(S/05,r) - expí- xt (1- € cos 8)] 
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Te cn the other hand P (S70; ) is fixed, 
O<P(S/0,r)<1, and the equation is solved for r, we will 
obtain an "egg-shaped" function with iso-levels Of 
probability of suppression P(S/0,r). The foxhole is located 
in the middle of the coordinate system. Along 9=0% the range 
r is a max for a certain fixed frobapility cf suppression , 
while at 9021805 r is a min for the same probability. . 


The runction takes the form: 


2 hn P(S/0 rT) 
rt = T € cos 8 (2.2) 


The follcwing graph shows the family of functions for 4 


representative selective prokabilities of sugffressions. 


E(S/@,r)=0.4°i i=1,2,3,4 (233) 


The parameters K and E are again assumed to be 1500 and 0.7 


respectively. 
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r/m = distance/meter 


| \ | 


r/n = o 
4 Z — 


foxhole 


nn I 
for 0-180 


ur 3 -  ISOBIBES 5° = =o Log 12 30. 
l-€ cos O 








— E 


These "egg-shaped"-functions of iso-probabilities of 
Suppression simulate the reaction of human beings lcoking 
along the main axis 0-0'in a very simple way. The functions 
take into account the visual and accoustical perception 
resulting from any given detonation of a round, where  ncise 
and light are the major stimuli. | 

One could think of other functions wich simulate tae 
behavior of an antitank gunner exposed te artillery rounds 


for instance: 


P(S/@,r) = L-exp{- = sl - e cos 08)] (2.4) 


This tyre cr function allows the probability of suppression 
E(S/0,r) to te smaller than one at its maxiual value and has 
the same general behavior as the functicn before. Asa 


enira modification: 
1 
P(S/8,r) s exp[- i dO - € cos 9)] (2.5) 


these functicns have the disadvantage that they have a 
discontinuity at r=0. The integration which is necessary in 
the follcwing derivation is more difficult than the chosen 
one. 

However, the crosscut secticns of these functions 
are not "egg-shaped" iso-functions but rather simple conic 
sections (ellipses) where the foxhole is located in the 
center of one focus point. 

The paper will continue with the function first 
described in (2.1), because of its simplicity and variety 
cf applicaticn. 

In order to derive a rate of suppression for the 


model, it is necessary to evaluate now, in a second step, 
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the unconditicnal probability of suppressicn P(S) for any 
incoming round, no matter where it will impact around the 
foxhole. This will be performed by matching the conditional 
probability P(S/0,r) with the area-hit-probatility P (A). 
Ccnsequently, to the engagement procedures of the 
artillery, the targets which shall be suppressed by tne 
artillery are categorized as small personal targets. This 
implies that the mean point of impact (MPI) of a given set 
of rounds lies on the target, which maans also, that the 
density has its max value at that point. According to the 
U.S. Army Field Manual FM 6~-161-1 page 2-2 it can be 
ascertained that the MPI-error is destributed normally with 
its mean at the alm-point. Thus it will be hypothesised that 
in such a case the distribution of inccming rounds is 


bivariate ncrmal with parameters 


= 0 (2.6) 


In crdar tc simplify the model, it is prcvosed that the 
dispersion of rounds expressed in tne standard deviaticn Is 


egual in all directions. With 


9 =0 (a7) 


ne correlaticn is assumed between horizontal and vertical 
deviation. Ihe normality is also preserved if more than one 
artillery gun is shooting. The essential change which has 
to be made when a whole artillery unit will deliver the 
rcunds will be the value of the standard deviation Ø 


In addition, o is determined by factcrs like: 
dispersion 


the fact that the location of the target is only 
estimated and an artillery unit delivers rounds in a 


fire area. 
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type of ammunition 


conditions o f the weapon systems (precision) 


ballistic properties of the rounds 

distance to target 

caliker p 
wind and other eee conditions. 


It follows frcm the bivariate normal density for artillery 
hits. 

— — 
a (Seer 


a 


£(x,>x,) = (2.8) 











2 (1-0°) 1 1 > 


and by inserting the previous mentioned assunntions and 


Changing to polar coordinates: 


2 2 2 
r = X1 + Xo 
x (2.9) 
1 
cos 8 = —- 
r 


that the unccnditional probability of suppression P(S) can 


be writteo in the form: 


AT rae 
J P(s/0,r)*£(r,0) de r dr (2.10) 
0 r=0 


0 
P(S) = 


|| Il 


0 


Recall that an area element in polar coordinates can be 


expressed: 


2 
r-dó-dr - d(5-) do 





Ihe conditicnal probability of suppression is 
P(S/0,r) = exp[- £ r" G - cos 6)] (2 11) 
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and the density for artillery hits f(r,9) is 





f(r,8) =—+5 expl- 5 e (2.12) 
270 ° 


inserting both into the expression above , we Can perfora 


the integraticn: 


= 2 
ES = f f expl- x sa - € cos 8) >, exp[- $5 Id) dé 
=0 r=0 210 (2,13) 





After a change of integration variable; 

















y? 
2 -=0 = 0 E RE 
2A 
2 I ( ) 
dC) - dz Aro 
we find 
08=2mw z=% > 1 
P(S) = f f exp [-($ a - e cos +). a d8 
8=0 z=0 g x. 
(2715) 
0=2T 1 
P(S) = J 2/2 1 
9=0 -270 E (1 - £ cos o) +7] 
g 
i^ Z =° 
en | - (2 (1 — £ cos o) ++} Jae 
K 2 
G Z=G 
o 1 
P(S) = d8 
8=0 het - E cos 9 + =z) 
2 Kc 
0-2 
K 1 
A z - de (2.16) 
410. 9=0 Ë +3) + (-e) cos 0 
2c | 
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Ey replacing 


Be + Ber (2.17) 
20 


the integral reduces to a known form which can be solved. 














K f do DIR 
P(S) = j where b >c 
2 ° 
E b t c*cos 0 
0-7 (b-c) tan 2 77 
K de K 2 2 
O pis ¿Ji A Er 
20^ 0-0 2:9 | /.2 2 (mme 
0 
P(S) = E. 2 ==> P(S) = Tou, ES) 


N 
| 
a 
N 
o” 
N 
j 
O 
No 
ia 


De Ve 


Inserting back the expressions for b and c, one obtains: 


A ss: (2.19) 
zt expe 
2c 


After the data analysis in chagter III B. and C. 
where the parameters K, € and oO will be determined, it 


will be possible to evaluate the probability of suppression 


29 





E(S) according to the preceeding formula. 

This probability adequately models those cases where 
temporary suppression is the only possible response tc a 
given detonation of a round. It may be useful to expand the 
model by introducing a second paired outcome at any given 


impacting rcund. Until now, we had one paired outcome: 


Ihe individual was either suppressed 
cr not suppressed. 
In additicn tc these we consider a second paired outcome: 


The individual is permanently suppressed by being 


wounded or killed 
or not permanently suppressed. 


If we focus the attention on modelling the erfects under 
this expansion, it is possible to evaluate a new probalility 
cf suppression P(S) and a probability of kill P(K) by making 


the following assumptions: 


The conditional probability of kill unlike the 
ccndtional probability of suppression does not 


deperd on tae aspect angle 0. 


It can be represented by a smooth curve of the 


Eclicving fora: 


P(K/r) = exp(- = ES (2.20) 


where H isa positive constant i.e. H> 0. This functicn was 
selected on intuitive grounds rather then based on 
real-world data. Taking a vertical cut thrcuch the surface 
function above along any angle f, one obtains the following 
figure. The parameter H is arbitrarily chosen as being 
H=50. Clearly this parameter is a function of different 


factors Jake terrain and amgunition. 
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r/m = distance/meter 


P(K/r) 


0.5 


r/m ? N ? i r/m 
for 9,2180? foxhole for 8,=0° 


1 2 
Figure-4 - FUNCTION P(K/r) = exp(- at ) 
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If on the other hand P(K/r) is fixed, O<P(K/r) <1, and the 
eguation is solved for r, one will obtain ccncentric circles 
with  iso-levels of probability of kill P(K/r) around the 
foxhole. 

Ine function has the form: 


* = - H:£n P(K/r) au 2.21) 


Ihe following Fig 5 shows the family of functions for "H 


representative selective probabilities cf suppressions. 
pia =0. 11 125225053. (2.22) 


Ihe parameter E is again chosen to be 50. 
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r/n 


for 0«0* 


10 
r/m 10 r/m 
8-90? 9=270° 
foxhole 


r/m 
for 0=180° 
r/m = distance/meter 
Figure 5 -  ISOLINES * = -H:2n P(K/r) 
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The concentricity around the foxhcle expresses the 
fact that no matter in which direction the person in the 
hole is looking the kill effect is just determined ty the 
distance. It would be beyond the paper to verify this 
assumption , and it would be extremely difficult to collect 
data for it. The reader must be content for the moment with 
the earlier presented intuitive argument. a 

This newly introduced function leads to a revision 
ct the prebability of suppression P(S) and the evaluation of 
the unconditional probability of kill, P(K). The density 


assumed earlier for artillery hits is used again: 





£(8,r) = —~y exp E: e QUIS) 


210 


Probability cr suppression Em) but not kill- is: 


8=27 r= 
PAG)- f. [f P@/e,r)[l - P(K/r)]-£(0,r) r dr de (2.23) 
0-0 r=0 
Remark: The subscript K at p, is used to indicate that this 
suppressicn probability appears together with tae 


PEODaADILity cf kill D(K). 
The expression P(S/6,r)-(1-P(K/r)) means that a round 
suppressed but did not kilil the individual at (1,0). 


Erobability of kill P(K) 


0=27 r>o 
PK) =f  f  P(K/r)-£(8,r) r dr de (2.24) 
0=0  r=0 


The computation for AS) and P(K) runs along similer lines 
as in the earlier evaluation of P(S). Fcr the specific 


functicns we find: 
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0=2T r=° 2 2 
1 f 
Ep et) * e exp[S O TOA) ae (2.25) 
020 r=0 210 
0227 rze 
2 2 
P, (S)= | f E r (i -ec cos 0)] [1 - m mM r )] 
K K H 
9-0 r=0 
2 2 
1 
> exp[- > © Id) do (2.26) 
2710 


change again the varlakle of integration 


” 





2 2 
r r 
1 - = Ü: zr 2-14 
5 z; CE 2.0 0; dS ) = dz (2.14) 
9=2 zo 
f 
P(K) = J f I 2 exp |- E + ix dz dé (2.27) 
6=0 z=0  2mo G 
E di 
P(K) = ——— (2.28) 
(G /H) + 1) 
9227 zo (2 1 1 
po) f f exp -[— (1 - e cos 8) + —]z — dz dé (2.29) 
K J ok 297 2 
0-20 2=0 G 2TG 


1 2 2 
— exp ¡7 (1 — £ cos 0) + dz de (2.30) 
G 


lI 
i 
No 


Recognizing that the first double integral is exactly the 


same as before we can simplify to 


8=27 I 
K 1 
x 2 i ER lee Be tane " 
K 2 Seg ee 3 
20 l + — - E G 
20 (2.31) 


35 











K f : de (2.32) 


P,(S) = P(S) - 
2 
K ATG 8=0 br -E Eje co cos 8 


The integration can be done by using the same formula as for 


the computation of PS). except that this time b is: 





b= (1+5 +8] (2.33) 
20 H 
From this it follows: 
2 
2 K K 2 
2 E] 
g Ë + 202 + 3 € 
ok 

PK 7 73 (2.35) 


A eu DT l 
20 2 2 
K 2 
Kee ea se 
20 ZO H 


After having derived the desired probapilities as summarized 
on the following figure, we are able tc evaluate the 
ditferent rates of temporary and permanent suppression. 

By assuming a certain fire rate A, with which the 
artillery 1S firing in the area where tne anti-tank gunners 


are located, the different rates will have tne following 


forn. 


Iwo-event model (unsuppressed-suppressed) 
Rate of suppression A, =A, PCS) (2.36) 
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Four-event-mcdel (unsuppressed-suppressed-survived-killed) 


Rate of suppression A = Ag PLC) (2:37) 


>> 
tl 


Rate cf killing 


k re P(K) (2.38) 


One important rate must still be developed; it is the Tate 
cf rise from a suppressed state. We model this by stochastic 


rrocess methcds. 


3,7 





P(S) = —— — — — (2.19) 
2 
jos i +5, n- 
20 
P(K) = — (2.28) 
O 
a + 1 
P (S) = P(S) = K (2.34) 


Figure 6 - SUMMARY OF THe DERIVED FORMULAS 
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2. Rate Of Rise 


The model so far represents a situation in which 
there are three possitle outcomes. When any independent 
round impacts, the ¡person in the foxhole is assumed to be 
either 

Suppressed, - 

nct suppressed, or 

is killed. 


Naturally this is only true if we assum2 that he is always 
back up again when the next round impacts. This fact krings 
up the necessity of considering the process along the tine 
axis. 

Eecause o£ the complexity of a stochastic process 
which reflects the tehavior situation cf an individual 
exposed tc arriving artillery rounds, it s¢ems useful tc 
start with a very simple process, in which the time of 
suppression ty a particular rcund is considered to be fixed. 
In a next step, this constant response time is randomized 
over the numrer of people under consideration, i.e. each 
individual has his own randcm time which hcwever is assumed 


constant for the process itself. 


Finally, the response time of a single person ray be 
considered to be a random variable coming from a certain 


distributicn. 


Both apprceaches are Simple, and represent a first attempt to 

escribe the actual situation in different ways. Certainly 
the suppression time can also be considered as a function of 
both above mentioned processes or of the miss distance r and 
the aspect angle 8, or of a combination of all. But these 
dependencies would rather complicate the mathematical 
derivations and lead beyond this study. 


Both processes consider rounds that arrive in the 
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neighborhood of the foxhole in accordance with a Poisson 
Erocess having a rate Ag. The model consists of two process 
states. One is the "up'-state which means that the person in 
the hole is unsuppressed and able to act according to his 
mission. The other state is termed the "down"-state which 
results from an incomming round and a possitle reaction to 
it in the form of suppression. In this stats the perscn is 


not able to fullfill his mission, he is physically dcwn. 


The first process assumes that the person in the  foxhole 
changes fror the "down"-state to the "up"-state only if he 
recognizes a gap of at least T time units before the 
appearance of the next round. This means that the 
"down'-state nas a duration of exactly time T if no rounds 
arrive in the time interval (0,T). The time T in connection 
with such a process may ke called the critical gap. 
(Saver,p.481]. The following Ligure shows the Easic 


relationship: 


4Q 





State O "up!" 
State 1 "dowm" 


— SE; t>T Tom orbutcal Pap 


¢ State 
Impact of first Impact of second suppressing 
suppressing round (process starts all 
round over agein 
T 
1 Å E a Bü 
06 $ 





| ; t/sec 
t — —— 





SECOND CASE: ¿<A 


A State 
Inpact ofurırst Impact of second sunpressing 
suppressing round (process starts all 
round over agein) 
t 





> 
A t/sec 
A u 


Figure 7 - SCHEMATIC DESCRIPTION OF THE MODEL 
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For the time being T is assumed fixed. However it is  nmcre 
reasonable to descrite it by a random variable, since tne 
time gap T is determined by factors like learning, 
accustoming to, cr overcomming of, fear, stress etc. The 
assumpticn of a fixed T has to ke changed in a later step, 
where it will be defined as a random variable. 

Furthermore, the change from the "up"-state to the 
“down"=state is accomplished by arriving rounds to which the 
person reacts through suppression ("down"-state) with tne 
earlier computed probakility of suppression P(S). 

Now let t be equal to tna total time of being 
continually suppressed. It is possible tc define T in the 
fcllowing way: 

T if no suppressing rcund arrives in (0,71) 

T= tl+t' if the next suppressing round falls (2.39) 


before T. 


In this casse the random variable t' is an inderendent version 


Cf the random variable +. 


According to this definition it is possible to compute the 
expected value of T which represents the mean time spent in 
the "down'"-state, or the mean time of being continually 
suppressed. The symbol used for it will be E{T]. 

In order to get the arrival rate cf suppressing 
rcunds which contribute to the time of being suppressed, we 
have to multiply the fire rate Ag with the grobability of 
suppression, P(S), which was evaluated in chapter II.B.1. We 
cbtain 


^s ^g PO) (2.36) 


First we have to state the following two expectations: 
E[tlno round in T] =T 


(2.40) 
E[rlround in T] = t!' 
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where t'= time of second suppressing round. Removing the 
condition on t' leads to expected time of being suppressed 


under the influence or incomming rounds. 





AT T At AT 
Elr]= Tte " * f e ^ a dt'[t' c EIX]] *E[s] ll-e Y] (2.41) 
0 
AT is Í AT 
Elt] = Tee SAE ESN 5 (- t - 1) + E[t] (l -e ) 
À 
S 0 
AT, SIT EUNT 
Er] Te " *z[1- Q*X,De *]+E[t](l-e *) 
1 -A T 
E[1] 2 T * P be (1 +A De ] 
TE 
y Ag? 
E[t] = = (e è - 1) >) 
ds ; 
inserting A = A et PCS) we obtain 
A_P(S) T 
1 £ =) (2.43) 


Plt] = YPG "° 


Extending the idea for a rixed time gap by applygiag 
it to a group of people separately, we are able to 


rerormulate the process in the following way: 


Assume a group of perscns in which sach individual sticks to 
a certain but randor time gap T when responding to a 
suppressiv= round, and assume that this random variatle T 


kased on the group is distributed with a density t (t) 
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The subscript P indicates the origin of the density 
(people). WE May express T which is defined as before as 


follows: 


t if no suppressing round 
arrives in (0,t) or t>t 
t =< where t is the time of the ° (2.45) 
first suppressing round. 
tti 3f tst 


R3mark,t' has the same distribution as To. A similar 


derivation as before leads to the following result. 


Elt|T=t] =t f Ae dt! + f [c'e E(x|lT-t]]A e ° det ’ (2.46) 

t 0 = 

At 

aecall, the expression Ae is the density of incoming 
suppressing rcunds. 

A E t =À t! cA t 
E{t|T=t] = te +f tie © de’ + EfrlT=tl(1-e °) 

0 

| ~À t ~À t 
E[t|T-t] - + 2 = e $ +Elrli=ti(l-e *) (2.47) 

S S 


Sclving for E{ r| T=t], we find: 


1 A 
Elr |T=t] = x (e - 1) (2.48) 
S 


cr replacing 4, again ty A, = Agr PCs) 


qu 





1 À eP (S)t 
E[t |T=t] = "I (e = 1) (2.49) 


So far this value represents the expected duration spent in 
the "down state" given a particular time gap t for a certain 
individual. Clearly the conditional expectation has the 
same form as the unconditional expectation E[t] (see formula 
(2.43)) for a fixed time gap T. 

Removing the condition in equation (2.49) we find 
the expected duration time of suppression based on the 


considered population. 


o9 


E[t] = | Elrl|r=t]- f (t)dt (2.50) 
0 
P 1 A ¿P(S)t l 
E[t] = — (e - 1)-f (t)dt (2:515 
0 X P (S) P 
| PPS 
1 f J 
Elt] = "IGI ^ e £, (t)dt = 1205) (2.52) 


For an evaluation of this expectation, the density £t) has 


to be known. 


A further variation of this precess leads to the 
second approach. Here the time gap T for a certain 
individual varies randomly according to a particular 
distribution. This approach is limited to one individual, 


and will not ke extended to a group. 
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Suppose an individual reacts with a random time gap t, 


T-t, 1-122234. 


cn any incoming suprressing round i, and assume that this 


random variable T is distributed with 


T ^ £ CO 1527539 


Ihe subscript T indicates the origin c the density (time). 


Ihe duration is defined as in formula (2.45) before. The 


T 
conditional ¢xpectation can be written as: 


-À t t -À t' 
Elt|T,=t] = te * +f (t' + £E{t]) ae 5 
1 0 S 


dt' (2.54) 


This time. T i5 conditioned on the individual's rirst chosen 


time gap I. 


l -À t =A t 
E(t |T,=¢] = (l-e ° ) + Elr](1- e E ) (2.55) 
S 


Removing the condition in equation (2.55) we can express the 
expected duration tine o£ suppression based on the 


individual's time gap distributicn. 


o9 


Elr] =f E(clT, -t]-£, (cde (2.56) 
0 
= 1 At -AE 

E{t] = f —(l-e *)+Er}(l-e *)/] £,(t)de (2.57) 


0 s 
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Solving tor the expected value E[t] and inserting 





= š 2.36 
A À ç P (S) ( ) 
we find: 
oo ALP (S)t 
f a-e ) f (t)dt 
2 (2.58) 
E(t] = oo -À P (S)t 
X P (S) Ë _ Í (1 — e ) £f Coa | 


cr expressirg in Laplace transform with s as an argument; 


E(] ^ 5 — - 1 (2259) 


gP (5) £(s) 





A further evaluation of this value reguires the distritution 
E(t) of T. 
T ) 


ccmparing koth expectations, 


Expected duration time of suprressicn based cn the 


pEcpulation 


co A P(S)t 
e £5 (t)dt - 1 (2.52) 


0 


1 
O 


and expected duration time of suppression based on 


the individual's time gar distribution 


Ss — | te 
Elt] = Ic) = 1 (2.59) 
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ke observe that they are different. However if we suppose 
that both distributions f(t) and f(t) are concentrated at 
T=z (delta function) "it 1s ne to reduce both 
expressions to the very first derivation (2.48) where the 


the time gap T was fixed. 


The reciprocals of these expectations approxinrate the rate 
at which the foxhole occupant (e. g. member of a group of 
antitank gunners returns from the suppressed state into the 


unsupgressad state) returns back to continue his mission. 





va E EET (2.60) 
whera E( 1] represents any of the derived expectations. 
In summary the three rate coefficients develcped are: 
A = As" PÈR) (2:38) 
— Xe P, (S) (2,37) 
e A (2.60) 


If killing as an additional event is considered, the rate of 
rise 1s computed with the same formula (2.60) except that 
this time BS) is used instead of P(S). 


ES SUPPRESSION MODEL 


The model which will be developed in this section can be 


set schematically in the following framework: 
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FORCE 
FIRE 
DUIS E 





NE EKK MOVING ARTILLERY 
GUNNERS TANK UNIT 


Figure 8 - SCHEMATIC REPRESENTATION OF THE SITUATION TO BE 


MCDELED 
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This analysis is restricted to the effects cn the anti-tank 
gunners... in this case the X-forces. It is behavior of 
members of the X-force that was modeled previcisly. 

The approach to formulate the situation cf the anti-tank 
gunners is done by using the ideas of lanchester-type 
equations and their further development { Tayler,1978,p.20]. 

The differential equations representing the model are 
all deterministic in the sense that each of them will always 
yield the same output for a given set of input data. Even 
though combat between military forces is a complex randon 
process. These equations shed light on combat dynamics and 
may be useful in defense planning studies. 

The kasic idea is that artillery forces use  "area'-fire 
tc suppress or eliminate forces like anti-tank gunners. 


"Area" in this context means the fact that 


the artillery unit "knows" the area in which to 
shoot, but does not know the location of each 


anti-tank gunner, 


tae anti-tanok gunners are "invisible" to the 


artillery unit. 


It we further assume hcmogeneous forces or X, it is possible 
to set up differential equations which model the rate of 


change Or the X- forces: 


dX (t) 
Ze EI RE) O CO b 0 AE) (2.61) 
dX (t) 

T = a X (t) - À X, (t) (2.62) 


Bemark: Clearly à here is taken from the Y-event-model. The 
S 
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variables X(t) and X(t) are respectively the number of the 
X-forces which are either active in the foxhole (able to use 
guns and anti-tank-weapons), X(t), or suppressed , X,(t).- 

In this simplified structure, tne two equations describe 
the most essential factors of the assumed situation. They 
are mathematically approximate , because the solucion of the 
tuo differertial equations will furnish the number of 
gunners active on the battlefield and the nugber of gunners 
supppressed at any point in time. The derivation uses the 


Laplace transformation: 


ri) CREE (0 
br (XxX(t)) 


(2.63) 


x(s) 


we perform Laplace transformaticn upon (2.61) and (2.62) 


sx_(s) - x (0) = = A + My? x (s) * À x, 6s) (2.64) 
sx (s) = x (0) = dx, (s) = A, 5 (5) (2.65) 
solving these two equations for x (x) and x (t) and 


translating them pack to X (t) and x (t) gives the desired 


time dependent quantities: 


-X (0) 
a 


- A. + A. + s) x (s) + A x (s) (2.66) 


-X (0) À x (s) - A, + s) x (s) (2.67) 


-X ON * s] 5 —0 £A. s) c^s)x (s)+1 QM ts) x (s) 
a S a u u S 


K 


li 


-X (0) Qc A A xa _ A, + s)A,, x (s) 


-X Ma, + s) - X (005, - A Aga CS) - A. | IA + s) x (8) 


x (0013 -X (0)2 F sX (O) 
S u a 


29 LU I ——-Ly 2.68 
x (s) = TX TX +S Ü *s)- AA ( ) 
e k u S u 
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multiplying (2. 66) and (2.67) with A and 


A + À, f s) respectively we receive: 


X (0)A. t XQ. + A) + sX (0) 


A AS N a 2.69 
x, (9) AU +2, + s)Q +5) - AA ( ) 
S k u S u 


Suppose that at t=0 the number of active gunners (able to 


watch and shoot) is equal to X. and the number of gunners 


0 
being suppressed (down in the foxhole) is equal to 0, 
for i.e. t=0 


X (0) = X 
a 





° (2.70) 
X (90) = Ü 
s 
the equations for X (s) and X (S) can be rewritten: 
X. À X s 
x (s) = ue — —e— (.7 
s + sa. +), + X) FAA, S + sa. +A, + A) u 
ZA 
x (s) = > e (2.42) 
S + síA. + A + A? + AA, 
Using the Laplace correspondence: 
Y es _ „pt 1 
A ¿CA p S for a #b 
s (-a-b)s + ab 
(2:73) 
pie u = , > 
== E or a 
a | * + (-a-b)s + ab 
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where a'b = u 
and -a- b A +A +À (2.74) 
S k u 
SEA +r. +r )J+AVOA_ +A p SEM 
(2.75) 
1 E 2 
and b fO +A + À 2] a Jo. + D + AD? - AA. 
[he equation for X (t) and X (t) then are: 
a S 
Xo at bt 
X (t) = [Q + a)e - QA + b)e `] (2.76) 
a a-b u u 
Xo at bt 
3 = = Lall 
Xx) ==% Ale e `] ( ) 
Since As A and À, have tc pe always greater or equal to 
zero 
A 2 OÜ 
— 
2,78 
120 ( ) 
ATRO 
— 
The constants a and b are always negative and real numbers. 
a < 0 
(2.79) 
b < O 
This leads to the basic shape of the functicn 
X (t) = £(t) 
à (2.80) 
X (t) = f(t) 
S 


shown in the following figure. 
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Figure 9 - FUNCTIONS X (t) AND X (t) 
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In discussing 


properties can be seen. 


X (0) 


x (0) 
s 


X (e) 
a 

X (°) 
s 


LOL t =t 
max 














the functions X (t) and X (t) tne following 
Xo a:0 b-0 
z ° es Zwol 
Zu CO, + a)e (A + b)e `] Xo ( ) 
Xo a+ 0 b+ 0 
A fe -e e (2.82) 
a-b S 
Xo at DE 
lim — [O,, Fa EP b):e jJ = O (2.83) 
t > œ 
X 
lim = le ne 120 (2.84) 
tt > ° 


X (t) has no max in 


X (t) 
S 


X (t) = X 


0 


has 3 max at 


, 


4 
S 


t 
max 


t 
max 


(1.e.nax values for X 


a 


(t) and X (t)) 


(0,9), i.e. max 


_ Qn (a/b) 
b- a 


Xo 


a - b 





À 


3 


With 


at 
max 


< 
[e 


value is at t€ = 0 


(2,85) 


(2.86) 


bt 
— < X, (2.87) 


The foregoing model assumed a rate of killing in additicn to 


a suppressicn 


rate 


and 


a rising rate. 


it is possitle to 


Simplify this model by leaving out the third rate. It will 


Le shown in chapter IIL C. 


k 


A, = AQ POO 


that 
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this rate 


(2.38) 





for such a scenario may be very small in ccmparison to the 
other two rates. 

For this reason we can state that leaving cut the 
killing rate will not drastically oversimplify the rodel, 
yet it will simplify the ccmputaticnal procedure In 
obtaining the wanted dependency between time and the number 
of gunners active or suppressed on the field. | ° 

This leads to the following set-up of differential 


equations: 


dX (t) 

ms -A,X (t) + dX Ce) (2.88) 

dX (t) 

aaa OA CE) =. E) (2.89) 
dt Sa u s 


where the total sum of people either suppressed or active on 


the battle field is equal to a constant X (no killing) 
Š U 


i.e. X zX (t)*X (t) 
Q a S 


which enarles us to rewrite the equaticns: 





š (2.90) 
Eam “A X(t) + A (X. - X (t)) . 

dX (t) 
F A (X. - X (t)) - A, X (t) (2.91) 


Using again Laplace transformation as earlier; 


À X 
a a 


S 





sX (s) -— Xy (0) 


=: = 2.92 
X X (s) + A.X (s) ( ) 


NAO 


s 





sX (s) - X (0) -AQX(S - AQX. (S) (2.93) 
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l 


Xo X (0) 


` = ———— + ——ə>— . — —ə 
x (s) sls+ (1 +1 )) s+ A +2) 
S u S u 
s sls + (A. + À )] s+ (A +2) 
S u S u 


using the Laplace transform: 


Il aeai] ado 
Liz wre ) í sís + a) 


-at 1 
L te a 


for both equations we find X (t) and X (t) fcr tne simplified 


(2.94) 


nodel. 
me. Tate A. Bas to De der yed tnis time. fron the 
probability c£ suppression P(S) without xilliag included as 


a possible event. We recall; 


Hugs ere I A =A, t P(S) (2.36) 
2 
nod 1 + ES s —* 
20 
Hence 
1 
= — - - - + 
X (t) A Xo +2, (1 - exp[-Q. + A,)t]) + XO) expl. A Del 
(297) 
1 
X (t) = A Xo x +h, (1 - exp lOS + A DED + X (0) exp[-Q +A Dt] 
for t = O | m) 
x (0) = X,, x (0) = 0 (2.97) 
a 0 S 
Tne eguaticn for X(t) and X(t) then are; 
a S 
Xo 
X (t) MED EE [A * n exp[- (4, + A t1] (2.93) 
S u 
F 
, = Y l = no = - 2.99 
X (t) Xo X (1 exp HA, 4 A De ( ) 
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Bemark: For t>% the sum of both assyagtotic values of 


X(t) and X(t) is equal to X. 
a S 0 


The basic shape of these two functions X(t) and X(t) can be 
seen in the following Fig 10. Í : 

If we ccmpare this result to the result on figure 9 we 
can see that the addition of killing to the model has an 
influence on the shape of the functions. So is e. g. the 
value for x (t) and X (t) in the first model (figure 9) for 
large times approximately zero, while in this case here 
(£igure 10), the total number of people on the field (X(t) * 
x (t)]) is always constant, i. ë. both functicns approaches to 


a limit value not equal to zero. 
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Figure 10 - FUNCTIONS X (t) AND X (t) 
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III. FITTING EXPERIMENTAL DATA DESCRIBING SUPPRESSION 


— — — — — — — — ——— — — — A A —— — — — — — — — —— — — — — — —— — —— — — — 


A. EXPERIMENTAL CESIGN 


1. General Aspects 


The mathematical model developed in the chapter tefore 

shall now be supported by an experiment which was conducted 
by the US Army Combat Developments Experimentation Command 
(Dac), Fort Ocd, California, 

The experiment descrited here belongs to a series of 
similar experiments, all dealing with tne objective of 
ccllecting and analyzing data of the suppression process. 
The data of this particular part of the CDEC experiment were 
based on the relationship between suppressicn and distance. 
This analysis tries to make use of the data by analyzing the 
relationship among suppression, distance, and angle. In a 
further experiment conducted by CDEC, the objectiv2 was also 
to include the angle as an additional variable for the 


suppressicn erfect. 


2. Setup And Realisation 


— — — —_— oe Se A 4 A eee ee SP nm 


The data were taken from a part of the experiment which was 
executed at Fert Hunter Liggett, California. 

Four foxhole-bunkers were constructed which 
guaranteed the safety of tte players as well as reproduced 
the real scenario as close as possible. Their tops were 
below grcund levsl, and covered with several layers of wire 
məsh and steel plates. This provided overtead protection 
from fragmentation. Each bunker was equipped with a mirror 


system which allowed the player to view events in front of 
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his position. 

A pop-up silhouette was installed fcrward of each 
player position. The player was able tc control the 
Silhouette as well as the cover of the mirror system, i.e. 
when the cover was opened (allowing the player to view 
down-range), the silhouette was in exposed posture. He was 
asked to respond as if he would be in the position of the 
pop-up silnotette. 

Each bunker was connected to controll bunkers by 
communication and instrumentation wires and power bunkers 
for data recording and supply. In the forefield cf the 
bunkers, different types of projectiles cr equivalent 
charges (31 mm,105 mm,and 155 mm among others) were placed 
and randomly detonated one at a time with the time between 
detcnations randomly selected from three possibilities of 
ten, firteen, and twenty seconds. The figure on the next 
page shows the schematic setup of the rounds, the location 
of the bunkers and the angle of sight for each foxhole. Tt 
can be seen that the explosions were visible to soue players 
but not to others. Since each type of ammunition has a 
different lethal radius, it was necessary to have different 
range configurations for each type. The aspect angle and 
the niss distances from the foxholes to the different points 
cf detonation are summarized in appendix A for each 
ammunition separately. Since all of the projectiles used in 
this part of the experiment were statically detonated, it 
“as not possible to model the kinetic contribution to the 
terminal effects. In order to keep the £ragmentation pattern 
as close as possible to those of incoming rounds, the 


projectiles were supported with sandbags [Suppex,p.A-1Zz]}. 
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The players were divided into two four-man teams. 
For each trial, memkers of one team occupied individual 
foxholes and provided the performance data for that trial. 
The mission of the players was to track moving target tanks 
by operating the periscopes. This mission was interrupted by 
the players responses to detonations in case he was 
suppressed (change of his state to  "down"). Thé state 
change and the pericd of suppression were automatically 


recorded. 


Appendix A summarizes the data which are the basis for the 
succeeding data analysis and serve tc evaluate the 
parameters and rates for the model. 


E. DATA ANALYSIS 


1. Parameters Of The Model 


— — — å P eee ee qm A — m F ib s — — 


The data presented in the section above are separated into 


three ditferent sets: 


Set cne consists of data related to the conditional 
frequency of suppression P(S/O,r) and the time of 
suppression when rounds of caliber 81 um were fired. 
Specifically, i£ n trials were made under conditions 
(1,9), the nunter s, c£ suppression was tabulated. 
Then s/n is an estimate of P (5/7, 9). 


Sets two and three provide the same data except that 
they are related to 105 mmm rcunds and 155 an 


rcunds. 


The main part of this analysis is to make use of the data in 
crder to estimate the parameters of the model described in 


63 





chapter II, and to see how well the real world situation can 
te described by the model developed earlier. 

The model is fitted (numerical values of the 
parameters are determined) by the method of least squares. 
De Groot, 668 With the basic model being of the 


form: 


E (3.1) 
Y Bu + Bv; + Err, 


And the data are analyzed on the following two different 


assumpticns; the second is certainly the most realistic: 


Homoscedasticity of the data, i. e. the error 


Variance of Err is assumed to be constant. 


Heteroscedasticity of the data, l. e. we will assune 
that each error term Err is distributed with 
variance c? "where the variance is not constant over 
observations. Errors are also assumed tc pe 
independent. If convenient, error terms will be 


assured to be approximately normally distributed. 


These assumptions imply different regression 
methods. For the first assumpticn, the result of unweighted 
single step and iterative regression methods will be 
presented. On the assumpticn cf heteroscedasticity, two 
Special iterative methods will be used. When we allcw for 
heteroscedasticity, ordinary least-squares estimation flaces 
the same weight on the observations which have small error 
variances as on those with large error variances. By 
applying a weighting regression, it is possible to adjust 
for the heteroscedasticity. So the two anncunced methods for 
the second assumption are iteratively weighted least squares 


regression methods. 
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In order to prepare the formula of the model fcr the 
regression, we apply a log-transformation to the equation 
(2.1) in chapter II B.1. This is a convenience, for it 
transforms the problem to one of linear fitting. We obtain 


the following result: 


Pitos = exol L EUER E OS (2.1) 
£n P(S/8,r) » - = : Er rs cos 8 (392) 


Using the notation 


In P(5/9,r) 


E 
l! 


u = r^ (323) 
V -Ir^cos9 


the equation can be applied to each datapoint i and the 


equation can Ee rewritten as: 


£ 
= m = 2 
Y; u k ( ) 


Y, =B u., + Bav., + Err (SL) 
l 2i 


(2.5) 


< 
o 
D 
H 
(D 
ct 
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m 
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and the randcm variable 


= $n P (3.6) 
Y, = £n P(S/0,,r,) 
The "hat" on the letter P represents an estimate of the 
probability of suppression at (r,8). The fitted regression 


^ 


plane for this function has the form: 


y = Bu + B,v Ä (3.27) 


The log-transfcrmation of tne conditicnal frequency 
cf suppression P(S/O,r) causes difficulties because of 
experimental results, which lead to an observed frequency of 
zero Suppression. This fact influences the sstinate 
P(S/8,r), that is used. 

Having specified the modal it is possible to apply 
the follcwing different regression methods cn the assumption 


of homoscedasticity: 


In other words, suppose n, observations were nade at 
experimental conditions (r, 9) and of these S were successss. 


Then first ccnsider the estimate 


Yi = gn P(S/r,,8,) = ın(s,/n,) (3.8) 


a G, 
where S; 


Of exposures to suppression, under condition i. 


is the number of suppressions, and n,is the number 


consequently, the transformed response to conditions (r,8) 
Day be ln (9/n,) = - © formally, causing embarassment. We 
will subsequently suggest alternative ways of dealing with 
this protlem. 
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Method 1: (zero-probakilities are omitted) 


With this method , the datapoints which had a 
frequency of suppression P(S/r,@)=s,/n, =0 are 
deleted. Consequently the following data points of 
the original sample of appendix A as shown in’ figure 


Fig 12 were nct considered. 
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Figure 12 - DELETED DATA EOINTS 
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Method 2 and 3 (Clustering of data points) 


One or more zero data points i , i.e. P(S/6, r) - 0 
i21,2,...0 Will be clustered with one other data 
point j which has a nonzero conditicnal frequency of 
Suppression, i.e. P (S/9,r) #0 and which possesses 
the closest distance proud angle SES the zerc data 


pcint(s). 


Py clustering (m*1) data points together and taking 
the mean of the frequencies , 8 and r within each 
cluster, we are able tc replace the data points of 
the cluster sets. This procedure is inevitably 


scmewhat arbitrary. 


A 1 m+l | 
P(S/0,r) * —T ) P,(S/0,r) 
i-l 
m+1 123 i mtl 15 i 


The bars above the symbols represent the cluster 
average. With these averages two different 


replacement prccedures can be applied: 
Each cluster will be replaced by its cluster 


average. This reduces the data sample. 


Each cluster point will be replaced by its cluster 
average. Conseguently the number of data  rcints 


stays unchanged. 


Figure 13 lists the replaced data points of the samples of 


appendix A. 
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Method 4 (data transformation by Cox) 


This method is analogous to one suggested by D. R. 
Ccx ([Cox, p.33]. The derivation for the variable y, 
shcwn on these pages can be applied  correspondingly 


in the following manner: 


S, +a ; 
J == (3.10) 
^J. ne 


where the random number S; is the number SE 


suppressions (successes, given n, trials). 


The constant a, which represents an unkiasing 
adjustment, is derived in appendix B. It is taken to 
Le; 


á => (1 - P(S/e,r)) . (3.11) 


This is suggested ty an auxiliary analysis similar 


to that of Cox. 


Hence the transformation results in the formula: 


1 ^ 
S, += (1- P.(S/8,r)) 
y = 1 ( i 2 — (3.12) 


ny 


This formula enables one to include data points that 
involve zero observed freyuencies cf suppression; 
the embarassment of taking tne logarithm of zero is 


no Longer present. 
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Method 5 (unweighted iterative regression) 


AS an extention of the method discussed last, this 


method replaces the value of P (S/8,r) in the formula 


S +5 (1 P, (S/0,r)) 


Y: = — U — — (3713) 
1 


After each iteration by an estimate of the 
probability of suppression P (s/8,r), the number of 
iterations was determined by the appearance of 


convergence of the parameters K and ce. 


The two rfcllowing methcds are weignted regressions which are 


necessitated ky the earlier described netercscedasticity. 


4 


Method 6 (weighted iterative regression Var (Err)er') 


This method is suggested by the fact that the 
ckserved variability in residuals increases with r. 
[ Pindyck, p. 100] Since the error variance is not 
known, it is reasonable to assume the existence of a 
Simple relaticnship between the error variances 
Var (Err,) and the values of one of the explanatory 
variatles in the regression model. In this 
analysis, the distance from tne foxhole to the 
explosion was chosen as the important explanatory 


variatle. By using 
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á 
Ve 4 Da 3 * 14 
ar[Err,] ry ( ) 


the resulting regression formula is: 


- 


Yi u; Vi — 
e deg n ee OT 
(Var[Err ])1/2 1 (Var [Err, 1)? 2 (Var [Err, 1)? 

Err, 
+ 1/2 
(Var [Err,]) 


This equation can te reduced to the following 


regression function: 


y; 7 84 * 8v; (3.16) 


where 


E e e 
in T 2 Í 


i 

Ca Ear 
EI E 
i 


(3.17) 


Afterwards, the iteration procedure is equivalent to 
method 5. This methed tends to make the variance of 
tne residuals around the fitted line of more nearly 
ccnstant variance; estimates or the model parameters 


should be thereby improved. 


(weighted iterative regression) 
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Using the error variance derived in appendix C, 


which has the forn; 


1 - P (S/0,r) 


Var[£n £i] = (3.18) 
n, P (S/0,r) 
where 
S 
£ =} (3.19) 
i n, 


the regression formula will be modified as follcws: 


Yi E Yi 
wa tB 1/2 * 95 172 
(Var[£n Eun (Var[£n £,1) (Var[2n f] 
(3.20) 
where y is again: 
1 A 
S. +ç (1 — P (S/08;r)) : 
n ee (3.12) 
i n, 


as derived in method 4. The resulting equation has 


the form: a S 
S. m (1 - P.(S/8,r)) bcp P. (S/0,r) 102 
M í 2 i zu 
n, : 
——-. = rw i T (3.21) 
€ = nn) v, 
1 Qt Err! 
u. * P. (S/O, T) ‘Ca a E 
n, * P, (8/8, r) 


The iterative procedure in this methcd is performed 
by using the estimate of the conditicnal probability 
of suppression as input fcr each succeeding 
iteration. The number of iteraticns was determined 


as in method 5 and 6 ty the apparent convergence of 
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the estimators of the parameters K ande. 


In order tc perform the multiple linear regression, we 
switch over to matrix notaticn, where we can write the 


e 
, 


normal equations in the follcwing forn: 
^ T , 
Z eZ-B = Z e. Y (3.22) 


The matrix Z is the design matrix which consists of vector 
U, the square of the distances and V, the product of the 
cos@ and the square of the distances. The variables y and Ë 
represent a (sample sizel*l)and a (1x2) vectcr respectively, 
For the experiments analyzed, the entries of the design 
matrix were cttained by solving the normal equations fcr 8. 
We find: 


8 = az try (3.023) 


Remark: Capital letters used for matrix notation. 


The parameters & and & can be evaluated by using tas 


equations: 
i 
By 


(3.24) 


m> 
u 


B7K 


As a modification of the suppressicr function (2.1) 
in chapter Il. B. 1 the function (2.4) in the same chapter 
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was considered. 


P(S/0,r) » Leexp |- z a - E COS 2 (2.4) 


Ihe inclusion of the parameter L allows for the possitility 
that suppression may not cccur for shots that fall in very 


close proximity to the foxhole. 


Method 8 (sizilar to method 4 but using model above) 


in this method the regression procedure and the 


method Y were appiied to: 


T: = Bu, + BV + 8, + Err, (3,25) 
where 
€ +4 a- a) 
i 2 i 
LX EU CE 
i 
Yi i 
_ £4 (3.26) 
Vi = ri cos ĝ 
1 
B=- 
€ 
Bo =X 
As 
The rcrmal equations have tne same form as befcre: 
B = (zT.zy l.;T.y (3.23) 


76 





but Z, the design matrix consists this time of the 
vectors U,V, and a column vector consisting only of 
ones. The resulting 8 vector consists now of 3 
elements, from which it is possible to evaluate the 


different parameters K, €, and L. 


The calculation of the parameters for all methods was 
pertcrmed fcr each type of ammunition separately by aprglying 
the computer language APL. The projrams which where used 
for this regression are listed in appendix I. The appendixes 
D through G show the analysis of the residuals after the 


regressicn cf the original data cf appendix A. 


The analysis for each ammunition type was performed 
according to the same pattern, and includes the following 


steps. 


For each earlier develcped methcd (1 tarough 7): 
Determination cf BETA(1) and BETA(2) 
Determination of K 
Determination cf € 


Plot of the residuals as a function cf y as defined 
in each method earlier. The APL function which was 
used has the name SCAT and belongs to the library 
package OA 2 3660 (available at W. R. Church 
Ccuputer Center). 


Plot of the residuals of the regression as a single 
array with the function BOXPLOT or the same library 
package. The plot characterizes the quartiles, the 
interguartile distance, the median, data points 
inside and outside the 1 and 1.5 interquartile 


distance and outliers [McNeil, p.13 and 71,72]. 


Numerical values of the residuals. 
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Histograms for the residuals of the regression which 
show the relative frequency and statistical 


features. 


For method 8, which is an analysis based on different 
assumptions for the prcbability cf suppressicn P(S/8,r), the 
same analysis was performed, with the exception cf the 
histograms fcr residuals. 

The methods 1-7 describe the attempt, to master the 
problem of zero-probabilities in the original data set as 
well as the problem of a possible heteroscedastisity. 

For the later methods, in which rore appropriate 
statistical tools were used, the systematic structure of the 
residuals seems to disappear up to a certain point. The 
residuals concentrate themselves more and more symmetrically 
around their mean and median (compare particular the ` 
boxplots and tha histograms in the appendixes D, E, F, and 
B). 

The relative large remaining range of the residuals 
is determined by single outliers. The analysis showed 
differences for different kinds of ammunition. Among the 
three ammunitions considered, the analysis cf 81 nm showed 
tne smallest spread for the residuals, for almost all 
regressicn methods. 

Tne appendix H shows for the iterative regression 
methods 5,6, and 7 the develcgpıent of K and € for the 
different kinds of ammunition. In each method , a 
ccnvergence of the values K and € with increasing 
iteration can be observed. The starting value for K is in 
method 5 smaller and in method 6 and 7 larger than its 
corresponding value after convergeance is obtained. This is 
true for all types of ammunition. An equivalent observation 
can be made for the € -values. 

À graph for the iterating € for method 5, 155 mnm 
could not te plotted because of the very small change cf the 


€ -values along each iteration. 
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The third page of appendix H displays the numerical 
values for K and € produced in each method for the three 
different types of ammunition. 

The analysis of the parameters K and € showed the 
fcllowing result; 

For the iterative methods (5,0,7) the K and 
e-values approached with increasing iteraticn a limit value 
(see page one and two of appendix H); this occured after 
about the sixth iteration. The final values for these 
rethcás are also presented in the tableau of appendix H 
(third page). The plcts for K and € show different shapes 
fcr different methods and partly also for different types of 
ammuniticns . The approach to the final value may occur from 
a relative small or a relative high value. Among the three 
anmunition types considered, the analysis of 81 am showed 
tne smallest spread for the residuals, fcr almost all 
regressicn ıethcäs. The scale parameter K influences the 
probability of suppression P(S/8,r) as stated under formula 


(2. 1)7 zn chapter II.B.1. In the following manner: 


An increasing value of K leads to an increase of the 
probablility of suppression F(S/9,r) .It is reasonable to 
assume that tne suppressicn effect increases with the 
increase of che calirer. This behavior was confirmed by the 
data analysis; except for method one and three where an 
inversion could be opserved between the K-values for €1 mm, 
105 mm, and 155 mm and for 81 mm and 105 am respectively. 

This distortion results from the fact that in rethod 
cne all data points with probabilities of suppression equal 
to zero are disregarded and in method three the cluster 
procedure emphasizes the average values reduced by the 
clustering. We prefer the latter methods. | 

Contrary to these observations on the K- values, a 
general trend for the e - values can not be related to 
different methods or different kinds of ammunition. The 
most reliable value of K and € seems to be found by method 
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7. This method considers the different variabilities for 
different sets of distance and angle. 

In the course of this paper, the value of K and € 
derived in zethod 7 are taken as input for developing the 
rates for the model of secticn II.B. 

The third parameter L developed by method 8 is valid 
only up to a probability of suppression evaluation equal to 
1.0. Since for 21 nm and 155 mu L is larger than one, the 
function for the mocel should be decomposed. This 
alternative approach to model suppressicn is beyond the 
scope of this paper and will not be considered further at 
this point, although it is certainly a topic worthy of 
further research. 

As a possible further step, confidence intervals for 


K and € could be developed. 


2. Suppression Tine 


The suppression tine data, which were collected in 
seconds, represent the time for which an individual remains 
suppressed as a reacticn of a single round. During this 
tine, the individual was unable to carry out ais missicn; ne 
is in hiding in an effort to survive. Accerding to the 
setup of the experiment, it was possibl2 for the suppressee 
to react to detonations, which he was able tc observe and to 
hear Or to hear only (visual/fauditory and auditory 
perceptors). This is consistent with the mcãel design in 
Chapter II. 

This analysis is an attempt to compare the 
suppressicn interval gained from the experiment with a 
certain distribution whose parameters are estimates of the 
data. The GAMMA distribution was selected, Lecause the range 
Cf the random variable X, representing the time, is limited 
telcw by zero and ty +” above. 

The data of suppression interval for the calibers 81 


mm, 105 mm, and 155 mm are plotted in histcgrams in apfendix 
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Ihe estimates A and r for the Gamma distribution 
are derived from the expectation and the variance of the 
data. | 


hith these parameters, the Gamma density can te computed in 


the considered range, by the formula: 


A 


ae ^ ^ 
x r-l -Àx 
a < < œ 
r() * : eS 
£(x/A,r) = (3.28) 
0 otherwise 


for compariscn, the related density was ccnguted by 


f) innis iatera (3.29) 
within the same rangs, where the frequency in a particular 
interval and the interval itself is taken fron the 
histcgrans of appendix J, and the constant N is the Sample 
Size of the measured suppression intervals. A . 

The numerical values for the estimates A and r and 
the values fcr the sample size N, the interval Ax, and the 


range are displayed in the fcllcwing table. 












INTERVAL 


^ 


Caliber E RANGE 


š 


81 



















Ihe following Fig 14 compares the Gamma density (3.28) 
with the height of the related estimated prcbability cf the 
data (3.29). The sketches show that the Gamma density for 
each ammuniticn type respectively underestimates the related 
frequency of the data because of its long tails. In fact, 
any distribution which meets the reguirements above could be 
taken for a comparison, although the Gamma distribution 
seems to ke a particularly good choice for the 81 mm caliber 


ammunition. 
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Figure 14 - COMEARISON GAMMA DISTRIBUTION AND DATA 
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In crder tc intensify the comparison, as a next step 
the skewness and the kurtosis of the data and of the Gamma 
distribution with the estimates above will be contrasted. 

These characteristics express the symmetry of the 
distribution about its point of central tendency and the 
relative concentration of cases at the center and along the 
tails of the distribution. 

For the derivation of the skewness and the kurtosis 
for a Gamma distributicn, we need the first four moments 


about the origin: 





oo -1 
r 2 2 -Ax (Àx)” 
— * = — = » = pS À d 
uu E[x] - E[x `] I! x e T(r) x 
(3.30) 
By changing the variable of integration 
z = Xx 
(3.31) 
dz = À dx 
we find 2 
2 " z^ -z r-1 1 1 s 
Uo = Elx ] = f 723 e z T(r) A Y dz (3,32) 
0 A 
u t „ırz-1 a? EFD i 
ATT) QO 
kut the akove integral is equai to one, hence 
NEUE br. (3.33) 
2 y 2 


In a similar fashion, the third and the fourth moment can oe 


derived. i 
Uy 7 3 + 2)(r + 1)r 
: (3.34) 
y — (r + 3)(r + 2)(r + 1)r5 
4 ,* 
Converting these moments to moments about the mean by 
using binomial expansicn: 
c 3 
x n 3 2 3:35) 
Uy = uz = 3, E{x] + 2(E[x]) | ( 
Y, = u, - Au, Elx] + 6n, (E[x]) - 3(E[x]) 


8y 





we can compute the skewness and the kurtosis 


— (3.36) 


For a Gamma distrikution we get the follcwing formulas. 


à GEHT 4 Hr, y, : 1 
3 3 2 A PONES 
À ) ES ny ° 
en 
3 yr 
and by Similar operaticn: 
= 3r + 6 
E" (3.37) 
Using the above estimates for A and r we get: 
3r + 6 
un = = ILL (3.38) 


2 
er Q 
e oe 


The formulas for the skewness and kurtosis derived from the 


data itself are: 


N 
1 ) - 3 
= I) iX, =x) 
oo 4 
3 N 
1 = 92970 
G È (x, - x)") 
N 
1 y =.4 
ST (x, — x) (3.39) 
en N 
4 N 
I =.2.2 
ES) (x, - 0) 
N t=] i 


The resulting numerical values for the sxewness and kurtosis 
derived fcr tne Gamma distritution and from the data are 


displayed in the following table: 
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81 mm ECEE ECEE | 155 mm | | 155 mm | 
Skewness a 


1.356 1.423 1.242 
GAMMA 
Skewness - 2.156 1.032 2.103 
DATA 
Rurtosis a, 5.759 EI c 5.315 
GAMMA 


Kurtosis La 8.004 1.309 6.750 
DATA 


In discussing the numerical features we can see that the 


L 


skewness of the Gamma distribution and of the data itself 
shows an asymmetric right-skewness (positive values). 
Although the values differ considerably . So is.e. g. the 
skewness of the data for 155 mm almost twice as big as the 
skewness based on the Gamma distribution. 

The same can be said for the kurtosis. There is 
significant difference between their numerical values. 

This analysis supports the fact that the Gamma 
distribution can only be a rough fit to the data given, and 
as already earlier expressed, a fit of another distrikution 


probably would have been as successful as this one. 


C. VERIFICATION OF THE MODEL 


In order to supply general features as input factors for 
the decision process for military leadership, the purpose of 
this paragraph is to compute numerical values for the 
probabilities and rates developed in chapter II.B.1. and 2. 

For the computation of numerical values, it is necessary 


to make the following reasonable assumptions: 


It is assumed that an artillery unit consist of six 
weapons all either 105 mm or 155 mm, or that three 
81 mg launchers are combined to a mission unit. 
Hence the fire rate A, for such units will be 
concluded to be [Wiener,p.189,211,213]: 


81 mr unit with 50, 55, and 60 rcunds/min 


105 am unit with 28, 32, and 36 rounds/min 
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155 mm unit with 18,21, and 24 roundsysmin. 


The standard deviation CG for the density for 
artillery hits will be considered to be [FM 6-161-1, 
Pad: 


25,30,and 35 m for 81 mm unit 
20,25,and 30 m for 105 mm unit 
30,3£,and 40 m for 155 mm unit. 


The numbers are taken from the field manual 6-161-1 
and out ot working papers of CDEC. They are rcunded 


fcr convenience. 


The parameter H of the probability of kill is chcsen 
completely arbitrarily with 100 for 81 mm, 150 for 
105 mm, and 200 for 155 mm and has nothing to do 
With experimentally cbserved values fer the weapcns 


here in questicn. 


The parameter estimates K and e€ are taken from 


apperdix H (tableau) with the values: 






1360.498 
0.871 


2349.845 
07399 


1991.151 
0.769 





0» »» 
l! 





Tne followinc figure displays the conditional probability of 


suppression with tne akove estimators for the parameters. 


P(S/8,r) = exp | - à ei - € cos 2 (3.40) 
K 
for 0=0 i.e. along the main cirection of sight. 
LL 3 
P(S/8 = O°,r) = exp l-> r (1 - o | (3,41) 
K 
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Range of r: -100 200 
P: O 1 


81 mm 
K = 1360.498 
Eps = 0,871 


105 mm 
K = 1991.15] 
Eps = 0.769 


155 mm 
R = 2349.845 
Eps = 0.859 


sees... ...... ......... 


a 
a 
One “a < ab “G sas q s < s s w x “q Qh an AAA O qu Q — 


° 
. e 
— — amp p — — —— — — — — — — —— s s s Q aw ow — y a — 6 


P. 
O een r 
P 
misc a .................... ..................... E 
p ° 

————— usss F 


Figure 15 - FUNCTION P(S/@ = 0°,r) = exp[- = oe] 
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For the above selected values of the standard deviation gy 


the fire rate \ and the paramet=rs K, c , and H the 


£ e 
values cr the probabilities 


A 


E sje A (3.41) 
^ 2 À 
20° ( +35) ~- me 
20 
^ ^ K 
pois DS SOTA (3.42) 
^ ^ 2 u 
20° IE - 5 
H 
20 

^ 1 

G /H + 1 


and rates 
A (with no killing) = re P(S) = 


7 Fr B à (3.44) 
A (with killing) Ae PAS) 


A = » P(K) 

K 
are displayed on the Succeeding two figures number 16 and 
17. Numerical values for the rate of rise A, are also 
computed and displayed in figure 17. The collected data 
(Suppression intervals) were not distinguished among 
individuals as in the course of the foregoing analysis. 


Because of this fact the numerical values for A, are based 


cn the expectation of duration time with regard to 
population. 
1 X P (S)t 
= —L———— = 35 
Eft] |} e £,(t)dt - 1 ( ) 


It might be worthwhile to evaluate the rate of rise 
by applying the Gamma distributions for t (t) which were 
discussed in chapter III.B.2, and contrast them to the 
values of A, Gained by the data itself. 
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Hence solving equation (3.45) by having set 





f(t) = e Q) à (3.46) 
ps? " 7%) Ä 
we find: | 
1-P(S)t At 
ERL £ e r-1 


if A > À P(S) and by changing the variable cf integration: 


(À - A ¿P(S))t = z ‘a X P(S))dt = dz (3.48) 
we get: 
r co =z r-l 
A e “ z 1 
E[r] = ———Əq—— — EON 5. 
APCS = aPC) T X P (S) 
— — — 
="] 
r (3.49) 
J 
À eP (S) [A ~ d-P(S)] £ 


Ey using the earlier developed estimators, we receive the 
Mathematical expression for the rate cf rise based on a 
Gamma distribution for the suppression interval data. 


^r -1 


4 = | —— A - (3.50) 
d-P(S)[A ~ A,P(S)]* en 
Renmark: This derivation is only true for A> AP(S) which means 
that for applying this formula, the fire rate À, nay not 
exceed the value; 


À 


Ctherwise the integral (3-47) explodes towards infinite 


and the resulting rate of rise would pe 
M20 — (3.52) 
u 
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which is consistent within the set of assumptions. 
Instead of using the Gamma distribution the data 
themselves were applied in a second ster to calculate an 


estimate for El(T] and À, based on the number of 


A¿P(S)t 
à s) - 1 13.33) 


observations n. 


A 1 1 n 

Er] = — = (à ) e 
PD "el > 

The approximate value for the rate of rise is then 


1 
E(t] 





i - (3.54) 
If killing as an additional event is considered, the rate of 
rise is conputed with the same fcrmulas (3.49) and (3.53) 
except that this time BAS) is used instead of P(S). The 
ccnputaticn of the values for the formulas (3.50) and (3.54) 
was performed by APL and FORTRAN respectively. The programs 


are displayed on appendix K. 
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PROBABILITIES P(S), p (S), AND P(K) 


81 mm 

g 
P 25 m 30 m 35 m 
P(S) O. OHIIIKCUNŠIIOIƏ 0.495 0.430 


P (S) 0.503 0.446 0.394 
P(K) 0.137 0.100 0.075 


105 mm 
20 m 25 m | 30 m 


155 m 


P 30 m 


P (S) 0.573 


P (S) 0.480 
P (K) 0.181 0.140 0.111 


Figure 16 - NUMERICAL VALUES FOR P(S), P. (S), AND P(K) 
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Figure 17 - RATES POR THE MCLEL 
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Row (4) and (6) and row (5) and (7) of figure 17 display 
comparable values for the rate of rise. As it was stated 
earlier, the Gamma distribution underestimates the related 
frequency of the data and hence the comparable values for 
the rate cf rise A, in row (4) are larger than the values in 
row (6). The same is true for row (5) and (7). They differ 
in general Ey 10 ~ 12X. If we are willing to live with this 
fact, the rates of rise evaluated by the fitted Gamma 
distribution are a good approximation fcr the values 
computed by the data. | 

In order to verify the four Lanchester equations (2.76), 
(2.77), (2.98), and (2.99), presented in chapter II.C., it 
is necessary to specify particular rates given in the figure 
before. In case of a certain known compcsition of fire 
units, specific rates could be developed as inputs fcr the 


model equaticns. 


By this, the model equations receive their specific shape 
and scale, their general behavior remains the same, as can 


Le seen in figure 9 and 10. 
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IV. ISCUSSION 


Based on the research effort of this paper it was fcund 
that the phencmenon of suppression as defined earlier, is a 
pultidimensicnal prcblen. It is influenced by 
psychological, physiological, and environmental variables. 
There exist many possibilities to model dependencies in 
qeneral form, however to make quantitative statements about 
suppression the modellers have to restrict their efforts to 
those variables which are observable and measurable. Since 
the main objective of this paper has been to establish 
models which are able to express relationships 
quantitatively, the main thrust has been to fcrıulate 
Suppressicn as functions of weapon systems and their 
dispersions. 

For the evaluations of the models, a set of simplifying 
assumpticns was necessary zn order tc guarantee a 
Dathematical transparency. Ihe dependencies developed in 
this tnesis postulate some satisfying results in modelling 


suppressicn. The models reflect sufficient accuracy of: 


The physionomy of the human being and its resulting 
behavior with regard to suppression. Suppression is 


mainly caused Ey visual and auditory perception. 


The influence of the weapon systzms i.e. their size 
and their firing capability are determining factors 


for the amount of suppression. 


Cf course the detailed results, i.e. the estimation of ncdel 
parameters, are based on selected weapcn systems and 
scenarios. Many possibilities exist fcr further work in 
this area particularly under the aspect of including 


extended human factors components in the construction of 
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suppression models. 


The authors feel the paper may provide a contribution to 
future design of wargames and Simulations as well aS weapon 
systems. It also supports a more careful analysis cf the 


combat situation. 
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APPSNDIX A 


ORIGINAL DATA 81 mm 
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APPENDIX A 


ORIGINAL DATA 105 mm 
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APPENDIX À 


ORIGINAL DATA 155 mm 
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APPENDIX A 


SUPPRESSION INTERVALS IN SECONDS 
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APPENDIX B 


DERIVATION OF CONSTANT A 


The ccnstant a can be determined correspondingly to the 


following derivation: 


Modifying equation (3.10) in Cox, The Analysis Of Binary 
Data, p. 33, we can define a transform as: 
S. +a 
yao ee | = (1) 


1 n. 
1 


Starting with the original model 


£(0,,r,) 
pa P.(S/8,r) = e 


and using the log-transformation 


— gn p. (2) 


we subtract 


9i + a 
uc dcc c O 


ae 
and chocse the constant a such that 


Fe = 1 (3) 
Ely, y ¿) 0 
Since S~Einorial (np, > np, (1 = p.) 


“hich can be approximated by 


S - n,P(S/0,r) * /n, U 
a l E 


where U is a randon variable with 
E[U] 0 


Eu?) p,(1-p,) 


Inserting S in equation (3) 
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en 


Ely, us y,] = £n P; 


MS 
_ u — 
U +— ) 
vn, P, iPi 


Approximating this by DE series 


mem a B EU our a f 
— n.p. 2 [A n.p. 
E p 1 1 i Pi l 1 
3 
— += me vet 
/n. p. iPi 
l L 


mt is sufficient to consider the first two terms of tne 

















Taylor expansion since U and a are small in ccaparison ton. 


= E U EE y? _ la _ BE 
a a 
1 1 1 1 l 


This equation has to be zero according to (3) 





Hence 
MM 1 
terms cr crder E(U] 
Yn 
i 
2 
terms of crder +. E a 0 
n p. 
T 1 2p, 


l 
e 


1 
29 (l-p,) = 


1 
a => (1-p,) 


terms of higher order are neglected. 
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APPENDIX C 


DERIVATION OF VARIANCE 


The weighting factor Var(ln f) for methcd 7 in chapter 
III B is computed as fcllows: | 

Expressirg the conditional probability cf suppression 
E(S/9,rT) as a quotient of the numper ot successes S. and the 


number of trials n, leads to: 


y 


A 
i n, 
l 

where the expectation cË f is: 


Elf.] = P.(8/9,r) 
i i 


S 
ín f, = £n (=) 
1 n. 


1 


S: 
Vari tn £.] = Var ES 
i : n. 


È 


Since tae random number S 1S binomial with mean 
E(S) = nP(S/0,r) 


in] 


and variance 
Var(S) » nP(S/0,r): (1 - P(S/0,r)) 


it can be stated as a function of a random variable U with 
mean 


E(U) = O 


and variance 


Var(U) = 1 
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= ° < e — > ` = . 

S, n; P,(S/0,r) T ^n, P. (5/0, r) [1 P,(s/e,r)] U 
Inserting this in the equaticn above: 

7 

> : ike N I 

n,P,(5/0,r) + /n,P,(5/0,r) [1 P¿(S/0,r)].U 


= 


[1 = P, (S/9,r)] 
= 7 e Á Ó 
Var[£n £ Var | £n P, (s/8,r) 1 + nP (870,7) U 


i = P Gona) ] ) 
HN NE D 
n.Po57/0,r) 
1 1 


Var[£n f] = Var | £n 


Var [£n £i] e War eh P, (5/0,r) - £n ( + 


Since P(S/9,r) is constant we knew that 


Var(2n P(S/8,r)) = 0 


which leads to: 


l - P,(5/8,r) 
Var[2n f] = Var £n l + a 70 e U 


knowing that for small x's 
Snel a km x for X $$ 


we can state 


~~ 


= P (S/O,r) 
Var[2n £1] ~ Var “a, P, S/O) s. U 


T 
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APPENDIX D 


ANALYSIS OF RESIDUALS 81 mm 
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ANALYSIS OF THE RESIDUALS 105 mm 
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ANALYSIS OF THE RESIDUALS 105 mm 
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ANALYSIS OF RESIDUALS (METHCD 8) 
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ANALYSIS OF PARAMETERS K AND 
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API REGRESSION PROGRAMS 


NON-ITERATION METHODS 
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[8] 
[9] 
[10] 
[11] 
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ITZRATION 
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REMARK: SIG CHANGES IN ACCORANCE TO THE 


DIFFERENT ITERATIVE METHODS, 
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HISTOGRAMS OF SUPPRESSION INTERVALS 81 mm 
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APPENDIX J 
HISTOGRAMS OF SUPPRESSION INTERVALS 155 mm 
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APPENDIX K 


PROGRAMS FOR COMPUTING RATES OF RISE 


Program for formula (3.54) 


FORTRAN IV G LEVEL 21 MAIN 
0001 DIMENSION TIM(500),SLA(20), RETAU (20) 
0002 N=348 
0003 READ (5,10) (SLA(I),I=1,18) 
0004 10 FORMAT (F10.5) 

0005 READ(5,11) (TIM(I),I=1l,N) 
0006 11 FORMAT(F10.5) 

0007 DO 1 1I=1,18 

0008 SUM=0.0 

0009 DO 2 J=1,N 

0010 CzEXP ( (SLA (I) *TIM(J))) 
0011 SUM=SUM+C 

0012 2 CONTINUE 

0013 A=SUM/FLOAT (N) 

0014 B=A-1.0 

0015 ETAU=B/SLA (1) 

0016 RETAU(1)=1.0/ETAU 

0017 1 CONTINUE 

0018 : WRITE (6,12) (RETAU (1) ,1=1.18) 
0019 12 FORMAT (1X,3F10.5) 

0020 STOP 

0021 END 


Program for formula (3.50) 


V LAMU 
[1] LA+(LAMSx ((LAM=-LANS >R ): (LALI*R ) 
[2] LAU+1: ((1:LA ) (1:LAMS )) 

V 
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